Structure of Bi2Sr2CaCu208+2; supermodulation from ab initio calculations 
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We present results of density functional theory (DFT) calculation of the structural supermodu- 
lation in Bi2Sr2CaCu2 08+a: structure, and show that the supermodulation is indeed a spontaneous 
symmetry breaking of the nominal crystal symmetry, rather than a phenomenon driven by interstitial 
O dopants. The structure obtained is in excellent quantitative agreement with recent x-ray studies, 
and reproduces several qualitative aspects of scanning tunnelling microscopy (STM) experiments as 
well. The primary structural modulation affecting the Cu02 plane is found to be a buckling wave of 
tilted CuOs half-octahedra, with maximum tilt angle near the phase of the supermodulation where 
recent STM experiments have discovered an enhancement of the superconducting gap. We argue 
that the tilting of the half-octahedra and concommitant planar buckling are directly modulating the 
superconducting pair interaction. 

PACS numbers: 74.72.Hs,74.81.-g,74.62.Bf,74.20.-z 
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The Bi-based family of high temperature supercon- 
ducting cuprate materials is alone among the several fam- 
ilies of cuprates in possessing an incommensurate struc- 
tural supermodulation. It is not known whether this 
wavelike distortion of the ideal periodic crystal structure 
represents a spontaneous symmetry breaking originat- 
ing in a mismatch between the preferred bond lengths 
in the perovskite and rock salt layers of the crystal, or 
is driven by disordered interstitial oxygen atoms which 
dope the system. Interest in this phenomenon, until re- 
cently viewed as a curiosity, has been reawakened by 
scanning tunnelling microscopy (STM) experiments [ij 
which show a strong positive correlation of the super- 
modulation phase with the superconducting energy gap, 
a local measure of the strength of electron pairing. Theo- 
retical calculations suggest that the structural mod- 
ulation may couple to local electronic structure, which 
in turn influences the effective local pairing interaction. 
Knowledge of how the lattice distorts in the presence 
of the supermodulation could be used to calculate local 
electronic structure and electron-phonon couplings, and 
guide a theory of local pairing in the high-temperature 
superconductors. 

The structure of idealized Bi2Sr2CaCu208+K , with 
two Cu02 planes per unit cell and a critical temper- 
ature of Tc w 9QK is shown in Fig. [H The half-cell 
shown is connected to an identical half-cell on the top 
or bottom translated with respect to the one shown 
along the a- or &-axis by half a lattice constant. The 
supermodulation is now a density wave superposed on 
the structure of Fig. (TJ characterized by a wavevector 
qs M which corresponds to a wavelength of approximately 
Asm = 27r/qsAf ~ 26 A k, 4.8 unit cells. The wavevector 
is perpendicular to the wavefronts along the a-axis, i.e. at 
45° with respect to the Cu-0 bond direction. Thus su- 
perconducting Bi2Sr2CaCu2 0g+a; is an almost-periodic 
crystal, since, strictly speaking, no two unit cells are 
identical fi^. Such structures were characterized math- 



ematically by de Wolff ^5:], and are known to occur in 
other materials such as Na2C03 as well. The supermod- 
ulation in Bi2Sr2CaCu208+2; leads to displacements of 
the atomic positions of up to 0.6 A from their ideal crys- 
talline locatioiis, as determined by a series of early x-ray 
diffraction 0, d, S [l3| and scanning tunnelling mi- 
croscopy (STM) experiments [Hj. 

Historically, the Bi2Sr2CaCu208+2; supermodulation 
has been regarded as largely irrelevant to high- 
temperature superconductivity because the material has 
a critical temperature of about 90K, similar to other 
cuprates with two Cu02 layers per unit cell which do 
not display this effect. Thus this type of lattice dis- 
tortion does not appear to have an essential role in ei- 
ther enhancing or suppressing superconductivity glob- 
ally. A good deal of attention has been paid to the ef- 
fect of macroscopic distortions of unit cell geometry on 
the critical temperature when external pressure is ap- 



plied [1^, |13|, In this case, however, exact changes 



in atomic displacements with pressure are impossible to 
determine because the materials are so complicated, even 
when structures are truly crystalline. On the other hand, 
a recent scanning tunnelling spectroscopy (STS) mea- 
surement [l| has shown that the Bi2Sr2CaCu20g+j; su- 
permodulation is correlated locally with the energy gap, 
one of the most fundamental observables characterizing 
the superconducting state. By carefully imaging the posi- 
tions of Bi atoms imaged on the surface layer of a cleaved 
Bi2Sr2CaCu208+a: crystal, the authors were able to asso- 
ciate to each pixel in the topographic scan a "supermod- 
ulation phase" ipsM, defining the relative position of each 
Bi relative to the maximum in the incommensurate den- 
sity wave. They observed a ^ 10% sinusoidal variation 
of the energy gap, determined by the half of the local en- 
ergy separation of the "coherence peaks" at positive and 
negative tip bias, with maxima generally coinciding with 
supermodulation phase 4>sm — 0. 

This measurement raises the intriguing possibility of 
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a new approach to understanding the origins of high- 
temperature superconductivity: by identifying in detail 
the atomic displacements which enhance superconduc- 
tivity, one may be able to determine the changes in local 
electronic structure or electron-phonon coupling which 
modulate the pairing interaction itself. While atomic- 
scale information is conventionally thought to be irrele- 
vant to superconductivity, these materials have coherence 
lengths approaching the nominal unit cell size, and there 
is evidence from STM studies to suggest that the pairing 
interaction can indeed vary on a scale of a few A [l^ in 
the presence of interstitial oxygen atoms. 

The key missing step in this program is the precise 
identification of the actual positions of the atoms present 
in the supermodulation of the complex BSCCO ma- 
terial. In principle, this information can be obtained 
by x-ray analysis, and considerable effort has been ex- 
pended [la . m, [131 in this endeavor. The analysis 
is hampered by the difficulty of analyzing superpositions 
of incommensurate harmonics [lit , as well as the presence 
of considerable disorder in the sample, thought to consist 
primarily of Bi/Sr substitutions on the lattice sites, as 
well as interstitial O disorder 2l| . The recent STM mea- 
surements do provide high-resolution imaging of atomic 
positions, but they are limited to the Bi atoms on the 
cleaved surface. 

In this work we report the results of density functional 
theory (DFT) calculations for the atomic positions in the 
Bi2Sr2CaCu208+a; system including the possibility of a 
supermodulation. While application of DFT to cuprates 
has been somewhat neglected because of its early failure 
to predict the insulating state of the parent compounds, 
it should produce accurate results for electronic states 
far from the Fermi level, and its success in making struc- 
tural predictions in a variety of complex materials is well 
known. 

Method. Details of the application of the DFT method 
to local structure problems in Bi9 Sr2CaCu20s+.T have 
been given earlier in and [23l |. The supermodula- 
tion to be studied in this work represents a wave of small 
deviations from the positions represented in Fig. [TJ In 
DFT treatment of solids, a small set of atoms is period- 
ically repeated. Such an approach is clearly impossible 
in the case of an incommensurate supermodulation, so 
we exploit the fact that the wavelength Asm of 4.8 nom- 
inal Bi-Bi unit cell constants of 5.4A is very close to 5 
cells. We therefore perform DFT calculations with peri- 
odically repeated set of 5x1 or 5x2 nominal unit cells, 
each of which involves 150 atoms (1420 electrons) or 300 
atoms (2840 electrons), respectively. By alternating the 
Bi2Sr2CaCu208+2; unit cell with an insulating blocking 
layer in the c direction, we can also study the BiO sur- 
face for comparison to, e.g., STM experiments. However, 
these calculations were performed for nominal half-unit 
cells only, due to computer time constraints. 

In addition, we study similar systems where a single 
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FIG. 1: Left: upper half of a nominal Bi-2212 unit cell. Lat- 
tice constants shown are between cations in each layer, values 
taken from Levin et al. [TH] (from this work). Right: side 
view of cation distribution for a layer of Bi2Sr2CaCu20g+i: 
one unit cell thick for (top) DFT bulk calculation, and (bot- 
tom) x-ray diffraction experiment [l^ . 



dopant O atom has been added. The lowest energy po- 
sition for the dopant in the nominal unit cell was deter- 
mined in Ref. [23], and shown to improve qualitatively 
the fit to the ARPES Fermi surface [111 . In the case 
of a supermodulation, the dopant may find new low en- 
ergy minima in each of the five inequivalent unit cells, and 
this calculation is therefore performed independently for 
starting positions of the dopant in each of them, allowing 
the system to relax in each case. 

Bulk results. We begin by noting that it is far from 
obvious that such a calculation will find that the sys- 
tem, in the absence of disorder, spontaneously deforms 
from its nominal structure shown in Fig. [T] to create a su- 
permodulation. To our knowledge this is the first DFT 
calculation which shows effects of this type (other than 
the much simpler ID Peierls-type distortions). Since the 
discovery of the supermodulation in Bi2Sr2CaCu208+a; , 
it has been controversial whether the phenomenon is in 
fact an unusual spontaneous deformation of a 3D periodic 
structure which lowers the total energy of the system, 
or whether the interstitials necessary to dope the super- 
conducting material stabilize the supermodulation. Our 
first result is therefore that, within DFT, the supermod- 
ulation indeed occurs spontaneously, as shown in Fig. [T] 
The waves of displacement from the nominal crystalline 
positions have significant amplitudes in all layers of the 
bulk system, and the largest displacements are in the xy 
plane, as also found in x-ray analyses. 

In the right panel of Fig. [1] we illustrate the remark- 
able success of these calculations by comparing a side 
view of the cation positions in the DFT calculation to 
those reported in the x-ray measurements of Levin et 
al. [l3|. The structures are nearly indistinguishable to 
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FIG. 2: Displacements in the BiO layer projected into the 
x~y plane for a) DFT bulk calculation; b) x-ray 16;]. Dashed 
lines are mirror planes, and X's mark points of stability found 
for O dopants by DFT. Dark green circles mark positions 
deduced by x-ray expt. c) Average z-displacements of Bi for 
DFT (solid) and x-ray (dashed). d)-f): same as a)-c), but for 
SrO layer. All arrow lengths were multiplied by 2 for clarity. 



the eye. Such a representation obscures the fact that 
the displacements in the x — y plane from nominal crys- 
talline positions are quite large and more complex, how- 
ever. These are displayed in Fig. [51 where they are indi- 
cated by arrows relative to the nominal sites. In the BiO 
layer, bulk DFT calculations and x-ray data are again in 
nearly perfect agreement. On the other hand, while cer- 
tain regions of the SrO layer show quite good agreement 
between experiment and theory, significant discrepancies 
appear elsewhere. This is because the atomic positions 
given by the bulk DFT calculations shown neglect the O 
dopant atoms, which were in fact located by Levin et al. 



at the sites shown in the lower panel of Fig. [2l The ef- 
fect of these dopant interstitials on the crystal structure 
was examined from the point of view of DFT both with- 
out [2^ and with the supermodulation present. The local 
displacements including the supermodulation are roughly 
as described in Ref. |23], but the supermodulation is found 
to break the symmetry of the effective potential for the 
dopant among the five nominal unit cells, leading to two 
effective potential minima 0.2 eV below other minima in 
the same layer, degenerate to within O.OleV, indicated 
by crosses in the figure. Thus we expect interstitial O 
dopants to populate these two different symmetry sites 
(shown in Fig. ^ randomly with close to equal probabil- 
ity. The existence of these two stable dopant positions 
is in agreement with STM [ij, but the x-ray work [l^ 
identified only one. 

Once the existence of these dopants in the crystal is 
accounted for, it is seen in Fig. O that the largest dis- 
crepancies in the position of the 0(Sr) (apical oxygen) 
between theory and experiment arise simply because the 
dopants push the apical O's to the side. The x-ray data 
display the displacements of this type averaged over the 
random positions of all 0-dopants. The dopants do not 
infiuence the BiO layer as much because the spontaneous 
supermodulation displacements move the stoichiometric 
O's near the dopant sites away from those sites. 

Within the framework of the DFT calculation, all qual- 
itative information about the all-important Cu02 plane 
is already given in Fig. [51 This is because in all our cal- 
culations the CuOshalf-octahedra are always found to 
be extremely rigid [23| , and tilt as a unit when an apical 
oxygen is pushed to the side as found here. Thus sig- 
nificant buckling (0(0. 4A)) of the CUO2 plane when the 
half-octahedra tilt is also a characteristic of the super- 
modulation. Note that the x-ray calculations do not find 
the same rigidity of the half-octahedra; we believe this is 
because the measurement averages over both filled and 
unfilled dopant sites. 

Surface results. In the lowest energy convergent DFT 
result for the free surface ("cleaved" at BiO layer) we 
have obtained thus far, the supermodulation structure is 
qualitatively similar to that obtained in the bulk case, 
with some minor surface reconstruction. For example, 
the amplitude of the z-displacement of the Bi atoms is 
enhanced somewhat, to about 0.57A between minimum 
and maximum, as shown in Fig. [3l This pattern should 
now be compared to the Bi positions obtained from STM. 
The smoothed Bi z-displacements from Slezak et al. [H 
are displayed for comparison in the same Figure: while 
the agreement is not bad, it is clear from examination of 
the topographic images in this reference and other STM 
works that there are also some fundamental differences 
between the x — y displacements observed by the two 
probes, suggesting a more significant surface reconstruc- 
tion than we have been able to reproduce by DFT. The 
STM pattern is found to correspond to a compression 
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FIG. 3: a) a displacements of Bi atoms in DFT surface calcu- 
lation (solid) and STM measurement [J]; b) schematic mod- 
ulation of measured spectral gap vs. supermodulation phase 
angle (f>sM [H; c) schematic tilting of CuOs half-octahedra 
and its dependence on (jisM. 



served by STM. This conclusion is identical to one we 
reached earlier in the case of O dopant atoms. The fact 
that two independent perturbations of the crystal lead 
to a) similar distortions of the half-octahedra, and b) 
similar local enhancements of the superconducting gap, 
strongly suggests that the electronic structure changes 
caused by these tilts are directly related to the local pair 
interaction. The obvious next step is to downfold this 
information to obtain local variations in hopping matrix 
elements, superexchange constants, and electron-phonon 
couplings. This work is in progress. 
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and rarefaction of the Bi positions, or alternatively to a 
wave in the Bi-0 bond angles, with no evidence for the 
Bi-O-Bi chain formation discussed by Levin et al. [3, El. 
We have found DFT solutions of this general type, but 
they are higher in energy, possibly because the half-cell 
calculations we have performed do not constrain the su- 
permodulation sufficiently. Thus we believe a DFT de- 
scription of the surface structure is possible, but requires 
further work. The goal is to complete the logical con- 
nection between the structural supermodulation, the gap 
modulation observed by STM (Fig. [8)3) , and the wave of 
tilted half-octahedra shown schematically in Fig. [S]:) 

Conclusions. We have shown that DFT calculations 
with unit cell 5 nominal unit cells long describe quan- 
titatively the displacements of atoms in the bulk of the 
Bi2Sr2CaCu208+2; material with its well-known super- 
modulation. To our knowledge this is the first ab ini- 
tio calculation able to accurately predict a 3D sponta- 
neous structural deformation of this type, and it provides 
strong evidence that the supermodulation is a sponta- 
neous breaking of crystal symmetry rather than being 
driven by interstitial defects. Calculations for a free BiO 
surface agree qualitatively in some respects with STM 
measurements, but differ in significant ways which we do 
not understand at this time. All calculations indicate, 
however, that the crucial structural element modulated 
in the 2212 case is the CUO5 half-octahedra; these ob- 
jects deform rigidly and in such a way that the max- 
imum tilt angle corresponds to the maximum gap ob- 
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